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Integral representations for three distributions of internal electric fields near isolated defects in 
ionic crystals are given. These three distributions are the Gaussian distribution for electric fields asso- 
ciated with phonons, and two Holtsmark type distributions for electric fields associated with charged 
impurities when the Jahn-Teller effect is small and when it is large. Numerical values for the distribu- 
tions and for the averages of squared-dipole matrix elements over them are tabulated in the conclusion. 
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1. Introduction 

The author and A. M. Stoneham [l] 1 examined 
recently the effects of internal electric fields produced 
by lattice vibrations and charged impurities on the 
lifetimes of low-lying relaxed states of isolated defects. 
We computed in [1] by statistical methods the distri- 
butions of internal electric fields produced by phonons 
and charged impurities. These distributions are the 
Gaussian distribution for the phonons, a Holtsmark 
distribution for charged impurities when the Jahn- 
Teller effect is small, and a distribution related to the 
Holtsmark distribution when the Jahn-Teller effect 
is large. We then averaged the square of dipole-radia- 
tion matrix elements over these three distributions of 
internal electric fields. We did not present in [1] the 
numerical methods by which we evaluated the distri- 
butions of fields and the averages of squared-dipole 
matrix elements over the distributions to avoid dis- 
tracting the reader from the major physical concepts. 

Some researchers have expressed an interest in the 
derivations of the expansions appearing in the appendix 
to [1] and in the numerical methods used to evaluate 
the averaged matrix elements. The author now takes 
this opportunity to present the derivations and numer- 
ical methods used in [1]. Such mathematical pro- 
cedures have applications to future studies on the 
temperature dependence of internal electric fields 
near defects. 

Section 2 contains integral representations for the 
distributions of internal electric fields and series expan- 
sions of the distributions which are valid for large 
and small values of the fields. Evaluations of the square 
of dipole-radiation matrix elements averaged over these 
distributions appear in section 3. Finally, section 4 
contains tables which give some representative numer- 



ical values for the distributions and the averages of 
squared-dipole matrix elements. 

2. Distributions of Fields 

The authors of [ 1] describe three classes of internal 
electric fields. The first class arises from the lattice 
vibrations. The second and third classes arise from 
the static fields associated with a random distribution 
of charged impurities in the host lattice. Some exam- 
ples of such charged impurities are substitutional 
aliovalent ions, vacancies, and interstitials. 

2.1. Lattice Vibrations 

Bennett and Stoneham [1] compute the probability 
that the internal electric field associated with optical 
phonons has a value between E and E+dE. They 
find that the probability for a harmonic lattice is pro- 
portional to the Gaussian distribution 



G(y) = 77- 1 / 2 exp(-7 2 /4), 



(1) 



where y = (2£2)~ 1/2 E is dimensionless and where 

G(y) is normalized to unity,! G(y) dy=l. The 

quantity £2 is the square of a characteristic electric 
field for the lattice. The Gaussian function G(y) also 
has the following integral representations; 

2 f x 
G(y)= — I exp (- x 2 /y 2 ) cos x dx (2) 

2 f x 
G(y) =- I exp (-y 2 ) cos (yy) dy. (3) 



and 



1 Figures in brackets indicate the literature references at the end of this paper. 



The representation (3) is useful for evaluating averages 
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of squared-dipole radiation matrix elements. Because 
eq (1) is in a convenient form, it is not necessary to 
derive expansions which are valid for large y and for 
small y. This is not the case for the remaining two 
distributions, however. 

2.2. Static Fields 

A charged impurity can perturb lattice defects such 
as F centers by a number of different mechanisms. 
One such mechanism is the electric field produced by 
the impurity. The authors of [1] also compute the 
distribution in the magnitude of the electric field |E| 
and the distribution in the projection of E onto a 
specific direction. 

When the Jahn-Teller effect is small, they find that 
the distribution in the magnitude of |E| =E is propor- 
tional to the Holtsmark function H(fi) , 



H(p) = — J °° dx x sin x exp {-(*//3) 3/2 }, 
ttQ Jo 



(4) 



where fi = a~ 213 E is dimensionless, and where H(fi) is 

normalized to unity, I H(fi) dfi = 1. The quantity a 213 

Jo 
is a characteristic electric field associated with the 
charged impurities imbedded in a dielectric host. 
Representation (4) is the convenient one from which 
to generate an expansion valid for large ft. The repre- 
sentation, 

H(P) = ^ | dyy sin (fiy) exp (-y 3 ' 2 ), (5) 

TT Jo 



also the useful one from which to evaluate averages 
of squared-dipole radiation matrix elements. 

When (3 has intermediate values near unity, the 
author finds that it is necessary to evaluate numerically 
the integrals in eqs (4) to (7). It is convenient for him 
to integrate by Gaussian quadratures with Legendre 
polynomials. The Legendre polynomials require a 
finite interval of integration. The semi-infinite interval 
^ x < oo is approximated by the interval 10~ 5 fi^ x < 
6.34 f3 when eqs (2), (4), and (6) are evaluated numer- 
ically in the present work. Expansions valid for (3 ~ 1 
coverage slowly and are not competitive with the above 
numerical integration procedures when computation 
time is important. Likewise, asymptotic expansions 
which converge rapidly when /3 is very large or very 
small are more efficient than the above numerical in- 
tegration procedures for such extreme values of /3. 
The next two subsections contain summaries of 
the derivation for such series expansions. 

2.3. Expansions for Large (5 

In this subsection, expansions for the distributions 
H(p) and //'(/3) which are valid for p»l are de- 
rived. The representation (4) for H (j8) and the repre- 
sentation (6) for H'(fi) are convenient ones from 
which to begin. When a = /3 _3/2 , the distributions 
become 



and 



H(p) = (2/7rp)I(a), 



H'(p)=(ll7Tp)J(a), 



(8) 



(9) 



where 



is the convenient representation from which to gener- 
ate an expansion valid for small /3 and from which to 
evaluate averages of the dipole-radiation matrix 
element. 

When the Jahn-Teller effect is large, the distribu- 
tion in magnitude of the projection of E onto a specific 
direction is proportional to a distribution H' (/3) 
related to the Holtsmark distribution; namely, 



dx cos x exp {— (x//3) 3/2 }, 



(6) 



1. 



H'{fi) = W~ 

7T/3 Jo 

where //'(/3) has the normalization^ //'(/})<#} = 

Both functions H and H' are called Holtsmark func- 
tions. The function H' is intermediate between a 
Lorentzian and a Gaussian distribution. Representa- 
tion (6) is convenient for developing an expansion valid 
for large j8 and the following representation, 



ff'W 



7T Jo 



dy cos (/3y) exp (— y 312 ) , 



(7) 



readily leads to an expansion valid for small ft. It is 



I(a)= dx x sin x exp (— ax 312 ) , 



and 



j( 



Jo 



dx cos x exp (— ax 312 ) . 



(10) 



(11) 



Because the quantity a approaches zero in the limit 
as ft approaches infinity, Taylor series expansions 
about a = become useful representations for the 
functions 1(a) and J (a). 

Using the following limits 

x sin x — — I — Lim 4- Lim I ( — e~ px ), 
2i \ P -*-i+ P^i+ 1 \ dp J 

and 

cos x = — I Lim + Lim I e~ px , 

2 \ p-+-i+ p-n + ) 
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p 



and the Laplace transformed function 

f(p; a) = I dx exp (— px — ax 312 ), 
Jo 

we write eqs (10) and (11) in the forms 

/(a) =V- (~ Lim + Lim ) r/tp; <*K 

and 

7(a) =-( Lim + Lim)/(p; a). 



(6m+l) 



1 a 2m+l » (6m + 8) (6m + 6) 

2(2m+l)! ik 2 2 



(12) 



(13) 



(14) 



(6m + 4) 



2 (2m + 2)! 



(17) 



15 /2 V' 2 24 



10395 
384 



©V 



11/2 _|_ 



(18) 



i? 



Expanding f(p;&) in a Taylor's series about a = 
and noting the Laplace transform 

dx exp (- px) x s = — , 

Jo p' +1 

where r(s-hl) is the gamma function, 

T(s+1) = \ dz z s exp (-2), 
J o 

and where s+l>0, we have 



and 



77/3 l\2/ £o 



(6m + 3) (6m+l) 



r/3 

1 a 2m+ 



2 (2m- 



1)! ^o 



2 2 

(6m + 6) (6m + 4) 

2 2 

2 a 2 »>+ 2 



2 (2m + 2)!j 



(19) 



/(a)=^(-Um + Um)2 



»'<*U 



(-a)" IZn 



T(f^) 



m + r) 



_3 *_» , 3g-« , 945 fi-™ 

L4(2-n-)" 2 n 192(27r) 1 / 2 



2.4. Expansions for Small fi 



(20) 



,,<3n/2) 



and 



y(a) = ^( Lim + Lim) J i-^ 



a) w V 2 / 



In this subsection, expansions for the distributions 
//(/3) and //' (/3) which are valid for £ <^ 1 are derived. 

Substituting y = x 312 in eq (12) and representing exp 
(— py 213 ) by the Taylor's series expansion 



n(3n/2)+l 



exp(-py^/3)=2 



{-pfl 3 ) n 



n=o nl 



we obtain 



Inserting — i = exp (— m'/2) and i = exp (tt/2) into the 
above equations gives us 



/(„)._ j izst r pa+ 



n=0 n! 



)(^.) -.(!=). 



2 • C 



-p)» U 3/ 

ft ! a [2(n+ 1)]/3 



(21) 



(22) 



and 



Applying the operations which are indicated in 
(15) eqs (13) and (14) to the representation (22) for/(p;a) 
gives us 



The expansions for Z/(/8) and H'(fi) which are valid i 

when /3 > 1, then become 

V()8)= AfM 1/2 f (6m + 5) (6m + 3) 
ttjS l\2/ Jf 2 2 
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sin 



2 °° \ 2 

/(«)=- y — - — — 

3~ 2 (b-1)! 



(^MM) 



a [2(n+l)]/3 



\3 3/ 



... 2 - cos (im/2) \3 V 
J(a)=- ^ 



3 n=0 



a [2(n+l)]/3 



(23) 



(24) 



501-968 O - 73 - 3 



Hence, the expansions for small /3 are 



m = v ' 



A. 

I™ 31T 



["-f^) + ^r(M) + ...] ; 



(*£*)• 



(26) 
(27) 



or 
#03) 

and 

2 °° /32m 

m = v 7 

or 

«'<»~£( r (I)-f+gr(f) + ...). « 

The first two terms of eqs (18) and (26) for the dis- 
tribution H((3) appear also in [3], To the author's 
knowledge, the expansions given above for the dis- 
tribution H' (/3) have not been reported previously. 



3. Averages of Squared-Dipole Matrix 
Elements 

Averaging the square of the dipole-radiation matrix 
elements over the distributions G(y), H(fi), and 
H' (fi) provides us with estimates of the effects which 
electric fields have on the lifetimes of relaxed F-center 
states. Discussions of these effects appear in [1]. It 
has been shown therein that the averages of the 
squares of the dipole radiation matrix elements over 
the distributions G, //, and H f are proportional to the 
following dimensionless expressions: 



M\ [G; y]=-Fl± P {1 + (a/y) 2 }- 1 ' 2 G(a) da\ 

(29 

Ml [H; /3] = ^[l ± j°° {1 + (a/p)*}- 1 ' 2 H(a) da}, 



(29) 



and 



(30) 



I [H';fj] = ± [l±2J X {l+( a l(3)z}-H*H'(a) da]. 

(31) 

The ordering of the energy levels of the F center 
determines whether the plus or minus sign is used. 
The details are discussed in [1]. 

The integrals appearing in eqs (29) to (31) are evalu- 
ated by two methods. Both methods employ numerical 
integration by Gaussian quadratures [2]. The first 
method is more straightforward than the second 
method; but it requires in some cases more time for 
computation than the second method requires. In 



the first method, the distributions G(a), //(a), and 
H' (a) are computed explicitly. Values of G(a) for 
the region ^ a ^ oo are given directly by eq (1). 
Values of H(a) and H' (a) are given by the following 
prescription. They are obtained when 1 < a <, 6 
by applications of Gaussian quadratures to the respec- 
tive representations (4) and (6). The respective ex- 
pansions (18) and (20) for H(a) and H'(a) are used 
when a > 6; and the respective expansions (26) and 

(28) are used when a < 1. These values for the dis- 
tributions G(a) , //(a), and H'(a) then are inserted, 
into the integrands of eqs (29) to (31). Finally, the inte- 
grals are evaluated by Gaussian quadratures. 

The second method does not involve explicit evalua- 
tions of the distributions. Instead, the representations 
(3), (5), and (7) are inserted respectively into eqs 

(29) to (31), namely, 



Ml[G;y] = -±- T da f" dy {1 4- (ct/y) 2 } 

2 77 Jo Jo 



-1/2 



M 



exp (— y 2 ) cos (ay), (32) 
i[H;P]=l±- f° da r dy{l + (a/py}-^ 

2 IT Jo Jo 

exp ( — y 3/2 ) ay sin (ay), (33) 



and 



Ml[H';P]=-±- J* da T dy{l+(alp)2}- 

2 7T Jo Jo 



1/2 



exp ( — y 3/2 ) cos (ay). (34) 

The integrations over a then are performed by re- 
ferring to several equations in section 9.6 of [2]. 
The results are 

T {1+ (a/0) 2 }- 1 / 2 cos (ay)da = /3K (l3y), (35) 
Jo 



and 



( X {1+ (a//3) 2 }- 1/2 a sin (ay)da = p*Ki(Py), (36) 
Jo 

where /3y ^ and where Ko and K\ denote, respec- 
tively, the modified Bessel functions of zero and 
first order. Insertion of eqs (35) and (36) into eqs 
(32) to (34) yields the following expressions: 



M 2 
and 



*[G; y]=~±l fcdyexp {-f)K {yy), 
U#;/3]=4±- frfyexp (-y^y/^QSy), 

2 TT Jo 



(38) 
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Ml[H';f3] 



2 77 Jo 



dyexp {-y*l*)Ko{p y ). (39) 



The modified Bessel functions are evaluated accord- 
ing to their polynomial approximations given in section 
9.8 of [2]. The final step of the second method consists 
of the numerical integration of eqs (37) to (39) by 
Gaussian quadratures. 

Eash of these two methods has its advantages. The 
decision as to which method is most suitable depends 
upon the values of y and fi and upon the other types of 
calculations included in the computer program. 

4. Conclusions 

We conclude this paper by presenting some numer- 
ical results. Table 1 gives values of G(/3), //(/3), and 
H'(p) for ^ /8 ^ 10.0. Table 2 contains the values of 
M 2 for each of the three distribution functions. The 
respective values of M 2 _ follow from the relation 
M 2 _= \-M%. 

Table 1. The distribution functions G(/3), H(/3) and H'(fi) 

All quantities are dimensionless. 



TABLE 2. The square of the dipole-matrix elements averaged over 
the three distributions G, H, and H' 

All quantities are dimensionless. 






Ml[G;y = (3] 


M\[H; p] 


M%[H';p] 


0.01 


0.51 


0.50 


0.51 


.04 


.54 


.51 


.54 


.1 


.59 


.53 


.59 


.4 


.74 


.61 


.73 


1.0 


.85 


.72 


.84 


4.0 


.97 


.91 


.96 


10.0 


.99 


.97 


.99 



P 


G(P) 1 


H(fi) 


H'(P) 


0.0 


0.56 


0.00 


0.29 


.04 


.56 


.00 


.29 


.1 


.56 


.004 


.29 


.4 


.54 


.06 


.27 


1.0 


.44 


.27 


.20 


4.0 


.01 


.08 


.01 


10.0 


.00 


.01 


.00 
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